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The transport coefficients for the Hofstadter model of non-interacting electrons in the presence of a 
magnetic field and strong on-site disorder are computed using an efficient numerical implementation 
of the non-commutative Kubo formula. The conductivity o and resistivity p tensors are mapped 
as functions of Fermi energy Ep and temperature T, and an asymptotic analysis in the limit of low 
temperatures is performed. In particular, the critical behaviors of the model at the plateau-insulator 
transition is simulated. The results reproduce all the key experimental findings on the plateau- 
insulator transition in the Integer Quantum Hall Effect: 1) The graphs of pxx as function of Ep at 
different temperatures intersect each other at a single critical point; 2) All these graphs collapse into a 
single curve after a one-parameter re-scaling; 3) The scaling exponents are in good agreement with the 
existing theoretical predictions; 4) The flow of o with the temperature, plotted in the {(7xy,(Jxx) plane, 
obeys the semi-circle law; 5) At the critical point, Oxx = ^xy = 6) The Quantized Hall Insulator 
phase, characterized by o^x - o^y - and p^y = h/e^, is observed at low temperatures. 

PACS numbers: 72.25.-b, 72.10.Fk, 73.20.Jc, 73.43.-f 



The plateau-insulator transition (PIT) and the plateau- 
plateau transitions (PPT) in the Integer Quantum Hall 
Effect (IQHE) are probably the most studied quantum 
transitions in condensed matter physics. The critical 
behavior of the transport coefficients near and at the 
transition points has been mapped w^ith extraordinary 
experimental precision, and various aspects were found 
to be in amazing agreement w^ith the single-parameter 
scaling theory [1 . 2] and the qualitative theory of IQHE 
130. However, some aspects remained puzzling, or 
became puzzling when confronted with new theoretical 
data, as we shall discuss below. On the theoretical front, 
to our knowledge, a direct simulation of the experiments 
was still missing, decades after the first high quality ex- 
perimental data on PIT and PPT emerged. This is a 
testament to the extraordinary difficulty of simulating 
the finite-temperature transport characteristics of disor- 
dered systems under magnetic fields. In this Letter, first 
and foremost, we want to announce that the numeri- 
cal approach of Ref. [6|, based on the non-commutative 
Kubo formula |[3|9|, enabled us to accomplish just that, 
at temperatures low enough to observe and analyze the 
scaling behavior of the transport coefficients at PIT and 
PPT (we present results only for PIT). Secondly, we want 
to announce that the simulations of the PIT with a dis- 
ordered model of non-interacting electrons reproduces, 
qualitatively and quantitatively, every aspect of the ex- 
perimental findings. These are important statements 
because, firstly, it means we have an effective tool that 
can be used to explore the remaining puzzling aspect 
of PI and PP transitions and, secondly, it demonstrates 
that the Quantized Hall Insulator (QHI) phase [10] ex- 
ists in models of non-interacting electrons with on-site 
disorder, contrary to what has been previously assumed. 



The Quantum Hall Liquid (QHL) phases correspond- 
ing to different plateaus in IQHE and the QHI pase 
are separated by phase boundaries containing extended 
quantum states [7J. As such, the PIT and PPT are 
essentially localization-delocalization transitions. The 
scaling of the localization length: ^ ~ (£f - Ep)~^ 
as Ep approaches the mobility edge combined 
with the concept of the temperature-induced effective 
system-size, introduced by Thousless yjj, and with 
the single-parameter scaling hypothesis |[T1 |2l, lead to 
a temperature-scaling law 1121 : 

p(£f , T) = F{{Ef - f ^)(T/To)-^), K = p/2v, (1) 

at the localization-delocalization transition. Here, f is a 
system dependent function and p is the the dynamical ex- 
ponent for dissipation. The finite-size scaling exponent 
V is expected to be universal (z.e. sample and model in- 
dependent), but the universality at PIT and PPT in IQHE 
is still debated, especially after the recent contributions 
that shook the old status quo, as discussed below. 

The temperature-scaling of the transport coefficients 
at PPT was resolved in the early works of Wei et al 11131 - 
15] and Kawaji and Wakabayashi |16|, where the scaling 
exponent k appeared to be sample-independent. Subse- 
quent experimental studies tT7'-'23] showed the opposite, 
but the non-universality of k was attributed to the dif- 
ferent dissipation mechanisms, and later to the nature of 
disorder [24]. The works by Koch et al |25^,'26l resolved 
both K and v, and the latter was in good agreement with 
the existing theoretical predictions at that time (v ^ 2.35) 
1271 . Subsequent measurements found similar agree- 
ment and by 1997 [28J the idea of universality in IQHE 
seemed to prevail. Same universality seemed to be at 
work for PIT lilQll29ti33 l. The old status quo was again 
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FIG. 1. (Color online) The density of states of the model Eq.|4] 
as function of the magentic flux O (in units of Oq = h/e). The 
horizontal white line shows the range of energies and the O 
used in our simulations. The QHL phase with Oxy = 1 and the 
QHI phase with o^y = are marked on the figure. In panel 
(a), W = 0, in which case one can easily recognize the classical 
Hofstadter butterfly. In panel (b), W = 3 and here one can see 
the impurity states partially or completely filling the spectral 
gaps visible at W = 0. 



reassured by the fine measurements from 2009 1341 per- 
formed on alloys with strong and short-range disorder. 
However, in the same year, a large scale simulation I35l 
found V to be about 10% larger than the previously ac- 
cepted value. As documented in Ref. f36l, this value of 
V ^ 2.6 has been reproduced (up to small deviations) 
by several other recent simulations II37H4T1 . The large 
discrepancy between the experiment and theory has re- 
ignited the debate about the universality in IQHE. To 
reconcile the experiment and theory, simulations of the 
experimental conditions are needed, to better under- 
stand how accurate is the ''saturation method'' used in 
Ref. 1251 l26l l34l to independently measure k and v, and 
to see if this method, in general, underestimates the ex- 
ponent V. We hope to convince the reader that this can 
be solved using the non-commutative Kubo formula. 

The second puzzling aspect is that the scaling expo- 
nent K was consistently found to be much larger at PIT 
(k ^ 0.58) than at PPT (k ^ 0.42) (m |33l [42-45]. We 
should mentioned that k and v were never indepen- 
dently measure at PIT, so at this point we do not have a 
direct confirmation of the universality of v between PIT 
and PPT. The different k values were ascribed to varia- 
tions in p due to different dissipation mechanisms. Here 
we report a v at PIT that is very close to the universally 
accepted value, so we can confirm this supposition. The 
third puzzle is the PIT transition into the QHI phase 
flO. 30, 45ti47l. In this phase, even though the direct and 
Hall conductivities gradually approach zero as T ^ 0, 
the Hall resistivity remains quantized ath/e^. From the- 
oretical point of view, the QHI phase is still a mystery. 
For example, in the global phase diagram (GPD) of IQHE 
derived in Ref. [48 J , the QHL phases are surrounded by 




FIG. 2. (Color online) The simulated (a) Oxy and (b) Oxx, 
as functions of Ep at different temperatures. The dashed line 
in panel (a) represents Oxy at T = 0, computed via the non- 
commutative Chern number. 

an ordinary insulator exhibiting the classical behavior 
Pxy ~ B/nec (in agreement with some early experimen- 
tal reports 1.49.-511). In spite of several theoretical efforts 
ll32ll52H59l , a consistent conclusion has never emerged. 
Examining the formula: pxy = Oyx/iol^ + o^y), one sees 
that the quantization of pxy is contingent to the peculiar 
scaling Oxy ~ o^^ as T ^ 0, which raises the following 
questions: Is this scaling due to a special form of dis- 
order (like long-range correlated disorder)? Is it due 
to the electron-electron interaction? We now can settle 
these questions: Since the QHI phase appears in our 
model, which includes only on-site disorder and does 
not include the electron-electron interaction, the answer 
is "no" to both questions. The observation of the QHI 
phase in our simulations was possible only because we 
were able to converge the finite-temperature transport 
coefficients at temperatures low enough to properly rep- 
resent the limit T ^ 0. 

Without further ado, we present the details of the sim- 
ulation and the results. We apply the non-commutative 
Kubo formula [7-9|: 

(J,y(T,T) = -r(fc,H](l/T + XH)"'[X;,OfD(H)]), (2) 

where, T represents the trace over volume, H is the 
disordered Hamiltonian, x is the position operator, t is 
the relaxation time for dissipation, Xh is the Liouvillian, 
and OfD is the Fermi-Dirac distribution for a given T 
and £f . The non-commutative Kubo formula was the 
key to the complete theoretical solution of the IQHE 
[7]. While Eq. |2] makes sense only in the strict ther- 
modynamic limit, we have recently [6] derived a canon- 
ical finite-volume (with periodic boundary conditions) 
approximation of^{T, t), which converges exponentially 
fast to the thermodynamic limit, as the size L of the sys- 
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FIG. 3. (Color online) (a) The simulated pxx as function of Ep, 
at different temperatures, (b) pxx as function of temperature for 
various Ep values. The arrow indicates the plateau-insulator 
transition. 

tern is increased: 

t) - cj^l:\Z t)| < const X (3) 

The convergence rate C,{T, t) was found |6| to be directly 
proportional w^ith T and inverse proportional with t, 
a conclusion that holds true regardless of the localized 
or delocalized nature of the energy spectrum at Bp. In 
real physical systems, t ~ T~P, p being the dynamical 
exponent appearing in Eq. [T] As such, provided the 
temperature or the system-size are large enough, our al- 
gorithm generates well converged results even near and 
at the critical points. A detailed discussion of the for- 
malism and extensive convergence tests for the present 
model, at the same temperatures and system-sizes as in 
the present study, can be found in Ref. |6|. A similar 
study for Quantum spin-Hall Insulators can be found in 
Ref. [60J. 

We use the Hofstadter model IBTI : 

H = -Y, e-'"*^"^'">c:c^ + WY, ^nclc„, (4) 

{n,m) n 

in energy units that set the hopping integral to 1. In 
Eq.jij O is the magnetic flux per unit cell (in units h/e), 
n A m = nim2 - n2mi, creates one electron at site n 
and (n, m) implies nearest neighboring sites. The co's are 
independent random variables uniformly distributed in 
the range [-1/2, 1/2]. The simulations are performed on 
a finite-size 140x140 lattice (hence L = 140), with W = 3 
and O = 22/L ^ 0.157. The transport coefficients are 
self-averaging, hence even a single random configura- 
tion can generate converged results if the system-size is 
large enough |6j. For the lattice-size used in our study, 
an average over five random configurations led to a root 
mean square comparable to the markers used in our fig- 
ures. Fig. 1 summarizes the spectral properties of the 
model and shows the energy range used in the simula- 
tion, chosen to cross from the QHL phase Oxy = 1 into 
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FIG. 4. (Color online) The flow of a with the temperature. Each 
flow-line corresponds to a fixed Fermi energy. The temperature 
decreases in the direction pointed by the arrows. The fine line 

represents the semi-circle: o^^ + (oxy - Ij) = {2^) - 



the QHI phase Oxy = 0. For clarity. Fig. 1 includes the 
W = case too. 

Fig. 2 reports the simulated conductivity tensor as 
function of Fermi level, at temperatures and dissipations 
kT = II T = 0.01,0.02,0.03,0.04,0.06 and 0.08 (k = Boltz- 
mann constant). Note that we are fixing p = 1 by making 
these choices. We have also included a calculation of the 
non-commutative Chern number (averaged over 30 ran- 
dom configurations), which is equal to cr^y at T = Q. 
The non-commutative Chern number was computed us- 
ing the numerical methods detailed in [62. 63]. In panel 
(a) we see the Hall conductance transitioning from the 
quantized value of e^/h to zero, and the transition be- 
coming sharper as T ^ 0. To a high degree, all the 
curves in panel (a) intersect each other at one point, the 
critical point. The slope of the curves at the critical point 
obey the well known scaling with temperature. Panel (b) 
reports the direct conductivity. For most of the energy 
range, Oxx displays an insulating behavior, with Oxx de- 
creasing as T ^ 0, except for a small region around the 
critical point where Oxx increases as T ^ 0, signaling the 
presence of extended quantum states. The data suggests 
that this region subsides to the critical point and that the 
maximum value of Oxx saturates at ^ ^ as T ^ 0, in line 
with the accepted theoretical picture of PIT IT2I . 

The existence of a single critical point, as opposed to 
a line of critical points, is best revealed in Fig. 3, which 
reports pxx = (^xx/i^^lx + ^xy) function of Bp (panel a), 
and as function of kT (panel b). Here, one can see (as ex- 
pected) Pxx 00 in the QHI side, but pxx ^ in the QHL 
side. While the latter behavior is usually attributed to a 
metallic behavior, in the present context it is due to the 
fact that Oxy ^ (see the pxx formula above). These oppo- 
site behaviors of pxx on the two sides of PIT are a bonus, 
because, unlike the Oxx curves, the pxx curves intersect 
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FIG. 5. (Color online) The simulated pxy as function of Ep, 
at different temperatures. The curves at lower temperatures 
display quantized values well beyond the critical point, which 
is marked by the vertical dotted line. For convenience we also 
show the data for p^x- 
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FIG. 6. (Color online) The simulated pxx as function of Fermi 
energy (a) before and (b) after the horizontal axis was rescaled 

as: £f ^ £^ + {Ep - E^) (^)^ with E^ = -3.15, kTo = 0.08 and 
K = 0.2. 



each other, very much like Oxy curves do in Fig. 2(a). Fur- 
thermore, to a high degree, the curves intersect at a single 
critical point, exactly as the experimental data show^ed 
(see Fig. 1 in Refs. |10|, |31| and |33|). The coordinates 
of the critical point can be extracted with great precision 
from Fig. 3(b) because, when plotted as function of fcT, 
Pxx displays distinct behaviors below and above a cer- 
tain critical Fermi level £^ ^ -3.15, with the lines curving 
downwards and upwards, respectively, as T ^ 0. From 
the same figure, we see that the critical value pj^ at PIT 
is virtually equal to h/e^, and from Fig. 2(a) we estimate 
a^^y to be virtually equal to 1 1-, exactly as in experiments 

Fig. 4 reports the flow of a with the temperature, in the 
{(^xyr (^xx) plane. With a remarkable precision, the flow 

is seen to approach the semi-circle al^ + (^Oxy - 

, and then to flow towards either Oxy = or 1, as 
T ^ 0. The semicircle law follows from general theoret- 
ical arguments in a "two-phase modeP' [64J of PIT and 
PTT. Experiments showed that the law is only approxi- 
mate at PTT t27J but it is extremely precise at PIT [10]. 
The reason for the latter is the fact that the semicircle 
law is equivalent to the quantization of pxy at ^ at low 
temperatures (see the expression of pxy above). 

The semicircle law observed in our simulations al- 
ready suggests the quantization of pxy in our data. And 
indeed, in Fig. 5, which reports pxy as a function of £f 
at various temperatures, one can see the pxy curves flat- 
tening out and converging to the line p^y = ^ as T ^ 0, 
exactly as observed in experiments (see for example Fig. 
2 in Ref. [45]). For example, pxy at kT = 0.01 remains ex- 
tremely close to this quantized value well beyond into 
the insulating phase. Three simulated quantized values 
occur when axy\T=o = 0.0005, 0.004, and 0.01, which in- 



dicate that the simulation box is much larger than the 
localization length. This and the semicircle law demon- 
strate beyond any doubt that, in the Hofstadter model 
with on-site disorder, the phase on the left side of the 
PIT is the Quantized Hall Insulator phase observed in 
Ref. [10]. 

The last piece of our analysis is the scaling reported in 
Fig. 6. Here we considered only the lowest four temper- 
atures because the scaling occurs in the asymptotic limit 
T ^ 0. As one can see in panel (b), the four pxx curves 
collapse on top of each other after a single-parameter re- 
scaling (see Eq.[l]). The best overlap of the rescaled curves 
was obtained with the scaling exponent k = 0.2 + 0.01, 
which compares extremely well with the value k = 0.192 
obtained from the formula k = p/2v with the presently 
accepted theoretical value v = 2.6 and p = 1 like in our 
simulations. 

In conclusions, the non-commutative Kubo formula 
enabled us to converge the finite-temperature transport 
coefficients of the Hofstadter model with on-site disor- 
der, at temperatures low enough to observe the univer- 
sal critical behavior at the plateau-insulator transition. 
The simulations reproduced, qualitatively and quantita- 
tively, all the characteristics observed in the experiments 
on IQHE. 
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